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Abstract. - We study the superconductivity in 2D fermionic systems near antiferromagnetic 
instability, assuming that the pairing is mediated by spin fluctuations. This pairing involves 
fully incoherent fermions and diffusive spin excitations. We show that the competition between 
fermionic incoherence and strong pairing interaction yields the pairing instability temperature 
Ti n3 which increases and saturates as the magnetic correlation length £ — > oo. We argue that 
in this quantum-critical regime the pairing problem is qualitatively different from the BCS one. 



In this communication we analyse the pairing problem in 2D fermionic systems near an- 
tiferromagnetic instability. Our key goal is to investigate whether or not the closeness to 
antifcrromagnetism is in conflict with the magnetically mediated d— wave pairing. This prob- 
lem is rather peculiar as on one hand the d— wave pairing amplitude increases at approaching 
the AFM instability due to softening of spin fluctuations 0] , while on the other hand, strong 
spin- mediated interaction destroys fermionic coherence and therefore damages the ability 
of fermions to form Cooper pairs. 

We demonstrate that the competition between strong pairing interaction and the destruc- 
tion of fermionic coherence yields a pairing instability at a temperature Ti ns which increases 
and saturates when the magnetic correlation length £ — > oo. We show that under certain 
conditions, Ti ns is universal in the sense that it does not depend on the details of the elec- 
tronic dispersion at energies comparable to the fermionic bandwidth W, and is determined by 
fermions located in a narrow region near hot spots - the points at the Fermi surface separated 
by the antiferromagnetic momentum Q. We assume in this paper that the Fermi surface does 
contain hot spots. 

We believe that the results of our analysis may be applicable to both cuprates and heavy 
fcrmion materials. For high T c cuprates, our results may be useful for understanding of the 
pseudogap physics in the underdoped regime, where the data show that the temperature when 
the system first displays superconducting precursors saturates at the lowest dopings We 
conject that our Tj„ s may be the onset of the pseudogap behavior, while the actual supercon- 
ducting transition occurs at a smaller temperature. For heavy fermion materials, our result 
may help understand the close correlation between the appearance of the superconductivity 
and an antiferromagnetic instability S. 
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Fig. 1 - Diagrammatic representation for the pairing vertex. The solid and wavy lines are fermionic 
and spin fluctuation propagators, respectively. 

The point of departure for our analysis is the spin-fermion model which describes low- 
energy fermions interacting with their collective spin degrees of freedom. This model can be 
viewed as the low-energy version of the lattice, Hubbard-type models, and is given by 

7t: = ^v F (k-k i .)4 Q c k!Q + ^X(7 1 (q) S q S -q + .9 4+q.a °"«,/3 c k ,/5 • S_ q . (1) 

k,a q q,k,a,/3 

Here Ck. a and S q describe fermions and collective bosonic spin degrees of freedom, respectively, 
and g is the spin-fermion coupling constant. The three input parameters in the model are the 
Fermi velocity vf, spin-fermion coupling g (which near half-filling is of order Hubbard U), 
and the spin correlation length £ defined via a bare static spin susceptibility which is assumed 
to be peaked at the antiferromagnetic momentum Q, i.e., Xo(q) = Xo£ 2 /(l + (q — Q) 2 £ 2 )- 
The dynamical part of the spin susceptibility comes from the interaction with the low-energy 
fermions and therefore is not an input. 

The model of Eq. ([!]) yields a spin-mediated pairing interaction which singlet component 
is r(q, fl) = 3g 2 x(q, 0), where x(q, 0) is the fully renormalized dynamical spin susceptibility. 
Near antiferromagnetic instability, this interaction is attractive in the d x %_ y i channel [H. 

A convenient way to study whether the spin-fermion interaction gives rise to a pairing 
at some Ti ns is to analyze a linearized equation for the fully renormalized d— wave pairing 
vertex F with zero total momentum and frequency. This vertex generally depends on relative 
fermionic momentum k and frequency to, i.e F = Fk(u)). In the ladder approximation which 
accuracy we discuss below, the equation for Fk(u> m ) takes the form (see Fig[l]). 

F k (u m ) = Fj®(w m )-T g J |^fc'Kj G k ,{w' m ) G-^-uQ T{k-k>,u m -u' m ). (2) 

Here Gfc(w) is the fully renormalized normal state single-particle Green's function. At T = 
Tins-, this equation should have a nontrivial solution even when F^\uj m ) = 

To analyse Eq.(||) we need to know the fully renormalized single-particle Green's function 
Gfc(w) and the pairing interaction T(q, Q) in the normal state. In 2D, the dimensionless 
coupling constant for Eq. ([!]) is A = 3g/(47r?;i?£ _1 ), where g = g 2 Xo is the effective spin- 
fermion interaction [||. Obviously, near a magnetic instability A > 1, and a conventional 
perturbation expansion is inapplicable. It turns out, however, that one can resum perturbation 
series and obtain a self-consistent solution for both Gfe(w) and T(q, Q) |§, |). This solution 
becomes exact in the formal limit N — + oo where TV = 8 is the number of hot spots in the 
Brillouin zone. Two of us have checked Q that the corrections to the spin-fermion vertex g 
are small by 1/N and can be safely neglected. 

The key effect captured by the self-consistent solution is the appearance of the small scale 
u> s f — 9/(8ttN) g/X 2 cx |g] which separates the regions of a Fermi liquid behavior at 
w,T < uj s f and quantum-critical, non Fermi liquid behavior at u>,T > uj s f. Specifically, for 
electronic states near hot spots, k w fc^ s , 

x(q,Q m ) - X o£ 2 /(i + (q-Q) 2 e 2 + ia»IM/) 
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G fe 1 (cj m ) — iu} m Zk(uj rn ) — ek (3) 

where 



Zk(uJ m ) = H 2 / ; ( 4 ) 



1 



Here e* = Vfc( k - kfes) and |vf(fc + Q)| = |vjr(fc)| = Vf- At T = and fc = khs, Z(uj m ) = 
1 + 2 A/(l + A /T+ w m |/w s/ ). 

Analyzing Eq. (|) at fc = fc^ s , we find that at w,T < ui s f, x(q, O) w Xo(q)j and G _1 (oj) 
has a conventional Fermi liquid form G -1 ^) w w + isignaj(w 2 + n 2 T 2 )/ (4u? g /). On the other 
hand, at w, T > 

-i/ - fq-Q^ 2 

X (q, fi) oc w 

V 9o 

G -1 ww + U<kT\ + (%\u\ w) 1/2 /(T/M)) signw (5) 



where u> = 4A 2 w s / = 9g/(2%N), qo — g/(2irvF), and f(x) is a smooth function with /(0) = 1 
and f(x ^> 1) « — 1.52^/ix. We see that spin fluctuations behave as gapless diffusive modes 
and fermionic excitations are fully incoherent. This behavior is obviously a quantum-critical 
one. Observe in this regard that Co does not depend on the spin correlation length. The scale 
Co will play a central role in our further considerations. 

The fermionic propagator also contains a linear in T term which does depend on £. This 
term, however, comes from thermal spin fluctuations which contribute n = m term to the 
frequency sum in Eq. (|4|). We will see that these fluctuations act as static impurities and do 
not affect Ti„ s . 

We first discuss in detail the pairing problem when (j/vf^f <C 1, i.e., when qo <^ hp- We 
argue that in this case, the pairing is dominated by fermions near hot spots and is insensitive to 
the system behavior at energies comparable to the bandwidth. Indeed, substituting the single 
particle Green's function, the spin susceptibility into Eq. (|^) and estimating the momentum 
integral using a d— wave condition Fk(uj m ) = —Fk+Q(uj m ), we find that typical \Q — q\ and 
\k — kh s \ are of order i-e., are much smaller than fc^- 

We also checked that for typical momenta, Zk(to) and Fk(to) are weakly k— dependent 
and can be approximated by their values at a hot spot, Z(to) and F(u>), respectively. Under 
these conditions, the momentum integration can be performed exactly. The N — > oo limit is 
particularly simple as typical momenta transverse to the Fermi surface are by a factor 1/N 
smaller than typical momenta along the Fermi surface. In this situation, the momentum 
integration is factorized: the one over transverse momenta affects only the fermionic Green's 
functions, while the integration over momenta along the Fermi surface affects only the spin 
susceptibility. Performing the integration we obtain 

F M = + AvrT^ . , (6) 

\u n \Z{uj n ) ^uj s} + \uj m - uo n \ 

Notice that the consequences of taking the N — > oo limit are the same as of the Migdal 
theorem for phonon-mediated superconductors: one can (i) explicitly integrate over momen- 
tum in the gap equation, and (ii) neglect corrections to g and to ladder series. More precisely, 
the 1/N smallness of the vertex corrections appears each time when these corrections involve 
fermions with momenta separated by Q B. For the spin-fermion vertex, this is always the 
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case, hence vertex corrections are small by 1/N. The pairing vertex has a zero total momen- 
tum, and the ladder diagrams for this vertex, which give rise to Eq. (Q), do not contain 1/N. 
However, the corrections to ladder series from, e.g., crossed diagrams do involve fermions 
with momenta separated by Q, and are small by 1/N. From this perspective, our analy- 
sis of the spin-mediated pairing is quite similar to the Eliashberg analysis for conventional 
superconductors 0. 

We now analyse Eq. (^). First we show that classical, thermal spin fluctuations, which 
account for inTX term in Eq. (Q), do not affect T; rls . These fluctuations account for the 
scattering with zero energy transfer and therefore act in the same way as impurities. Ac- 
cordingly, our argumentation parallels the one which shows that nonmagnetic impurities 
do not affect T c in conventional superconductors Introducing F m = F(uj m )/r] m where 
rj m = 1 + (XttT '/ 'Z(io rn ) /\u> rn \), we explicitly rewrite Eq. @ as the equation for F m 

F m =F^+XnTT — Si ^ (7) 

where Z is the same as in Eq. (|J) but without the contribution from m = n term in the 
frequency sum. We see that Eq. (Q) contains only the contributions from quantum spin 
fluctuations. 

We next discuss the form of the kernel in the r.h.s. of Eq. (]?]). We see that it contains 
two energy scales: u s f cx and £— independent ui 3> w s /, which is the upper cutoff for the 
y^uj behavior of the fermionic propagator. For \uj\ > Q, the kernel converges as 1/u 3 / 2 , i.e., 
the pairing problem does not extend above u>, which for g < vpkp is still much smaller than 
the fermionic bandwidth. 

The presence of the two energies u) s f and uj raises the question on how Ti ns depends on 
£. To address this issue, consider the form of the kernel in Eq. (^) at different frequencies. 
At |a->| < ui s f, the system behaves as a Fermi liquid (Z(u>) « 1 + A). In this frequency range, 
the kernel reduces to a constant, i.e., the pairing problem is of BCS type, with the effective 
pairing coupling constant X/Z = A/(l + A) which never becomes large. If frequencies above 
uj s f were not contributing to pairing, Ti„ s would be of order ui s f e~' 1+A )/ A , i.e., it would scale 
with u> s f . This is similar to what McMillan obtained for conventional superconductors |J . 

Consider next \lo\ > uj s f. Here the pairing interaction (the last term in the r.h.s. of Eq. 
(0)) becomes frequency dependent and gradually decreases compared to its zero frequency 
value. At weak couplings, this decrease obviously makes frequencies larger than u s f ineffec- 
tive for pairing. However, at large A the situation is more tricky both in our case and for 
phonon superconductors ]To|| . The point is that for large A, the mere reduction of the pairing 
interaction above uj s f is not sufficient - one also has to neutralize the large overall A factor in 
the r.h.s. of Eq.(§). At u> < Lu s f, this overall A is neutralized by Z(uj m ) « I + A. However, 
above uj s f, Z(u) m ) decreases as Z(u> m ) ~ A(w s //|w TO |) 1 / 2 , and the effective coupling X/Z(u> m ) 
increases. Simple power counting shows that this increase exactly balances the decrease of the 
pairing interaction such that the 1/ |oj | form of the pairing kernel survives up to frequencies 
of order uj. This may sweep the pairing instability to a temperature Ti„ s ~ uj ~ g/N. 

To illustrate this point we introduce a dimensionless parameter ut — (u> / '{nT)) 1 / 2 and 
consider the limit u> s f 0. In this limit, Eq. ([?]) simplifies to 

p = + - V Pn " r (8) 

m m 2 ^ ^2|n-m|vl2^ + Tf n T + ^+T\ 

A fictitious parameter a (= I in our case) is introduced for the subsequent perturbative 
analysis of this equation. We see that at low temperatures, i.e., large n T , the kernel in Eq.(pl) 
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has a 1/n form typical for a pairing problem. On general grounds one might expect that the 
pairing instability occurs at tit = 0(1), i.e., at Tj ns ~ ui |nj|. If this is the case, then the 
pairing is dominated by frequencies where the fermionic excitations display a fully incoherent 
quantum-critical behavior, i.e., the pairing is qualitatively different from that in a Fermi liquid. 

The above argumentation is, however, only suggestive as it is a'priori unclear whether Eq. 
(^|) has a nontrivial solution for any tit- Indeed, on one hand, the l/u> n form of the kernel in 
Eq.(||) is typical for a pairing problem and gives rise to the logarithms in the ladder series. On 
the other hand, this kernel depends not only on the running frequency as would be the case 
for BCS superconductivity, but also on the frequency transferred by the interaction. This last 
frequency serves as a lower cutoff for the logarithmical behavior. 

To get further insight into the problem we assumed that Fm is a constant and analyzed 
Eq.(||) for various a. We found that for small a, when perturbative analysis of the logarith- 
mical series is valid, the dependence of the kernel on the transverse frequency is crucial, and 
even at T = 0, the summation of the series of logarithms give rise to a power- law behavior 
F m oc /\iu m \ a / 2 rather than to a divergence. In other words, unlike BCS theory, at a <C 1, 
the logarithmical series do not give rise to a pairing instability. 

We find, however, that the convergence of the perturbation theory is confined only to small 
a -c 1. Indeed, assume that at small w m , F m oc |cj m | _1 / 4+/3 . Substituting this into Eq.(||), 
we obtain an equation on (3: 1 = (a/2)$(/3), where 

_ 7T^2 1 1 

$( ' ~ y/2 r(3/4 + j9)r(3/4-/3) cos tt/3 - cos tt/4 ' (9) 

For real f3, is an even function of (3, which increases monotonically from $(0) w 8.97 

and diverges at [3 — * 1/4 as $(/?) ps 1/(1/4 - 0). For a <C 1, we find (3 = 1/4 - a/2, i.e., 
F m oc |cj m |~ a / 2 , in agreement with the results of the summation of the logarithmical series. 
As a increases, (3 becomes smaller and reaches zero at a = a cr — 2<I>~ 1 (0) w 0.22. At larger 
a, a solution with real (3 is impossible, i.e., a perturbation theory breaks down. Instead, the 
condition 1 = (a/2)<&(/3) yields an imaginary (3 = if3* i.e F m oc |w m | -1 / 4 cos {[3* log |w m |). 
Near a cr , we find f3* ~ 1.2(o — acr) 1 / 2 . The appearance of the oscillating solution at T = 
implies that the pairing susceptibility is negative for some \u) m \. This obviously signals that the 
normal state at T = is unstable against pairing. An estimate of T;„ s may be obtained from 
a requirement that a temperature should exceed a maximum frequency where the pairing 
susceptibility is negative. For sufficiently small (3* this yields Ti ns oc u) e -7 ^ . We see 
therefore that for a = 1, when (3* = 0(1), the attraction between fully incoherent fermions is 
capable to produce a pairing instability at Ti ns ~ Co ~ g/N, as we conjected above, but this 
result has a non perturbative origin. We also performed RG analysis of the leading 1 /N vertex 
corrections and found that they only slightly, by 0(1/ N), change a cr which still remains much 
smaller than 1. 

To check this analysis, we solved our original Eq. (Q) with F$ = numerically for various 
A. The results are presented in Fig. (||). In the limit A — > oo we found Ti ns w 0.17cD. It 
is interesting to observe that the weak dependence of Ti ns /Cd on A, which is an indicative 
of quantum critical superconductivity, persists down to A ~ 0.5. This means that even at 
moderate A the pairing instability has a non-Fermi-liquid, quantum-critical origin. 

We now discuss the momentum dependence of Fk(u m ) at 2j ns . This momentum depen- 
dence is likely to mimic that of a pairing gap at T < Ti ns jL2j. As we said above, F(u m ) 
along the Fermi surface is weakly k dependent at relative deviations from a hot spot by less 
than g/vpkF which is a small parameter in the theory. We checked that at larger deviations 
from a hot spot, F(uo m ) rapidly decreases, as l/(k — kh s ) 2 ■ This means that for quantum- 
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Fig. 2 - The results of the numerical solution of Eq ([?]) for different values of the coupling constant 
A. 

critical pairing, the d— wave pairing gap is more strongly confined to hot regions than a simple 
cos k x — cos k y form. This result is intuitively obvious as the very fact that the pairing problem 
is confined to hot spots implies that the pairing state is a superposition of many eigenfunctions 
from the B\ g representation with almost equal partial amplitudes. Simple manipulations with 
trigonometry show that in this situation, the slope of the gap near the nodes should be smaller 
than the one inferred from the gap value at hot points assuming cos k x — cos k y dependence of 
the gap. Notice, however, that this effect is non-critical, i.e., the width of the gap in A:— space 
remains finite even if £ = oo. 

Finally, we briefly discuss the situation at large spin-fermion interaction, when g 3> vf^f, 
i.e., qo 3> kp (see Eq. (^). In this limit, the momentum integration extends over the whole 
fermionic bandwidth, and the presence of hot spots at the Fermi surface becomes less relevant. 
The explicit evaluation of T ins is no longer possible, but the reasoning along the same lines as 
above shows that T ins is independent on £ and scales as the largest typical frequency for the 
pairing problem. This typical frequency is obtained from the condition that maximum \q — Q\ 
are of order kp, and is obviously J ~ («f&f) 2 /(Ng). 

The analysis of the system behavior below Ti ns requires one to solve a set of three coupled 
integral equations for the fermionic self-energy, the anomalous vertex, and the spin suscepti- 
bility. Setting this aside for a separate publication fl3[| , we merely argue here that the pairing 
state which emerges below Ti ns is highly unusual and has no analogs in BCS superconductors. 
Indeed, on one hand, Tj„ s and hence the gap at T = are independent on £, on the other 
hand, the resonance frequency of the spin mode scales as 0J res ~ ^f£ -1 ~ Ti ns / X and for 
A 3> 1 is much smaller than the pairing gap. In this situation, it is tempting to conject that 
superconducting coherence may be destroyed by fluctuations not included in the Eliashberg 
treatment at T c < Ti ns , yielding a disordered region between T c and Ti ns . This issue is, 
however, highly speculative and requires further study. 

We now briefly discuss the situation in cuprates. Near half- filling, wf^f scales with the 
fermionic bandwidth, while g is of order of the Hubbard U , hence J and Ti ns (if g 3> if&f) are 
of order of the exchange integral of the corresponding Heisenberg model. The actual situation 
in cuprates probably falls into an intermediate regime g > wf^f- We emphasize however that 
for u/ g f ~ 10 — 20meV, and A ~ 1 extracted from NMR experiments at optimal doping 
the universal result (Fig. @) yields Ti ns ~ 10 2 — 10 3 K which is a reasonable estimate. The 
non-critical sharpening of the superconducting gap with underdoping is also consistent with 
the recent photoemission data Jl5| . More detailed analysis requires a more precise knowledge 
of both A and u> s f for various doping concentrations. 
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Finally we discuss how our work is connected to earlier studies. The Eliashberg-type 
equations for magnetically mediated pairing have been analyzed several times in the liter- 
|, mostly using the numerical technique. In particular, Monthoux and Lon- 
zarich Wl\ recently solved Eliashberg equations for large £ and for the Fermi surface with hot 
spots. They found that for large couplings, Ti ns likely saturates at a finite value at £ = oo. 
This fully agrees with our result for Ti ns . However, our key finding is the discovery that near 
antiferromagnetic instability, the pairing problem is a quantum-critical one, and is qualita- 
tively different from the BCS pairing. We also found that in the presence of hot spots at the 
Fermi surface, T ins is universal and does not depend on the form of the pairing potential at 
lattice scales. This physics was not detected in earlier works ]l4 
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